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Abstract 

In these lecture notes I present a short review of nuclear shapes, shape coexistence and 
shape-phase transitions in the interacting boson model. In a study with random interactions 
it is shown that the appearance of regular spectral features is a far more common phenomenon 
than was previously thought. The origin of these features are explained by studying the 
relation with the underlying geometric shapes. 

1 Introduction 

The concept of shape coexistence and shape-phase transitions in atomic nuclei is currently 
a topic of active research both experimentally and theoretically ^ The phenomenon of 
shape coexistence refers to the occurrence of different shapes in a given nucleus and has been 
observed in fission isomers in the actinides, intruder orbitals in or near closed (sub)shell nuclei 
and superdeformed nuclei at high spins Whereas coexisting shapes are usually associated 
with different configurations {e.g. normal and intruder), shape-phase transitions describe how 
one geometric shape evolves into another within the same configuration, e.g. between spherical 
and deformed shapes in the Nd, Sm and Gd nuclei or between 7-soft and deformed in the Pt-Os 
mass region [1]. 

It is the aim of these lecture notes to present a short overview of shape-phase transitions in 
the interacting boson model (IBM) and its extensions and to discuss the appearance of regular 
spectral features from the IBM with random interactions and its relation with the underlying 
geometric shapes and critical points. In the first lecture (Section 2), I present a review of shape- 
phase transitions for the IBM including two more recent developments, i.e. two-neutron transfer 
reactions as a probe of shape-phase coexistence in the Ge isotopes and the shape-phase diagram 
for the neutron-proton IBM. In the second lecture, I discuss some surprising results that were 
obtained in studies of the IBM with random interactions: the predominance of ground states 
with angular momentum L = and the occurrence of both vibrational and rotational structures 
despite the random nature of the interactions. 

*Lecture notes: IV International Balkan School on Nuclear Physics, Bodrum, Turkey, September 22-29, 2004 
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2 Shape-phase transitions in nuclei 



Shape-phase transitions in atomic nuclei were studied extensively in the early 80's in the frame- 
work of the IBM. The general procedure was laid out by Gilmore in the 70's using atomic 
coherent states [S] in combination with catastrophe theory ^ which was applied to nuclear 
physics by Dieperink, Scholten and lachello PP, Feng, Gilmore and Deans [7] and Lopez-Moreno 
and Castafios |8j . After the introduction of critical point symmetries as special solutions to the 
Bohr Hamiltonian [^l^j, the subject of shape-phase transitions has witnessed a revival both 
experimentally and theoretically. For a recent review of shape-phase transitions and critical 
point symmetries in nuclear physics, see P]. 



2.1 The interacting boson model 

The IBM provides an elegant and powerful tool for the description of collective nuclei In 
this model, collective excitations in nuclei are described in terms of a system of N interacting 
monopole (s^) and quadrupole {dl^ with m = 0, ibl,ib2) bosons with = and 2"*", respec- 
tively. The the number of bosons N is determined by half the number of valence nucleons. In 
addition to N, the eigenfunctions have good angular momentum and parity L^. Let's consider 
the schematic IBM Hamiltonian of the consistent-Q formulation (CQF) [12) 

H = ehd-KQ{x)-Q{x) , (1) 
where fid is the number operator for quadrupole bosons 

hd = V5{d^ X , (2) 
and Q{x) denotes the quadrupole operator 

Qmix) = (s^ xd + d^x +X{d^x d)'i^ ■ (3) 

Here s = s and dm = {—l)'^~'^d-m- The Hamiltonian of Eq. describes the main features 
of collective nuclei: it contains the dynamical symmetries of the IBM for special choices of the 
coefficients e, k and X) and allows to describe the transitional regions between any of symmetry 
limits as well. 



2.2 Potential-energy surface 

The connection between the IBM, potential-energy surfaces, geometric shapes and phase tran- 
sitions can be investigated by introducing a coherent, or intrinsic, state which is expressed as a 
boson condensate [3] 

|Ar,/?,7> = -L=(&t)'^|0) , (4) 



with 



= ^7^= (st + /3cos7 4 + -^/3sin7(4 + dL2)) , (5) 
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where the variables (3 and 7 determine the geometry of the nuclear surface. Spherical shapes 
are characterized by /3 = and deformed ones by /3 > 0. The angle 7 allows one to distinguish 
between axially deformed nuclei, 7 = 0° for prolate and 7 = 60° for oblate deformation, and 
triaxial nuclei 0° < 7 < 60°. 

The potential-energy surface is given by the expectation value of the Hamiltonian in the 
coherent state jl] 



y(/3,7) = {N,P,^\:H :\N,P,^) 
eiV/32 kN{N - 1 



4/3' - 4x 



^/33cos37 + ^x'/3' 



(6) 



The normal ordered form of the Hamiltonian : H : was used to eliminate the contributions of 
the one-body terms of the quadrupole-quadrupole interaction to the energy surface V{(3,'y). 

The equilibrium shape associated to the IBM Hamiltonian can be obtained by determining 
the minimum of the energy surface with respect to the geometric variables P and 7, i.e. the first 
derivatives vanish 



dV 
'dp 



dV 



, 



(7) 



and the determinant of the Hessian matrix of the second derivatives is positive 



a/32 



dj2 



I 5/3^7 ] 



> 



(8) 



2.3 Dynamical symmetries 

First some general properties of the IBM are reviewed, in particular the equilibrium shapes for 
the dynamical symmetries that arise for special choices of the coefficients e, k and x in the 
Hamiltonian of Eq. A more detailed discussion can be found in jllj . 

The ?7(5) limit 

For K = 0, the Hamiltonian H of Eq. reduces to the ?7(5) limit of the IBM 

Hx = ehd . (9) 
The potential-energy surface of Hi is given by 

MP) = ■ (10) 

The equilibrium value of the deformation parameter /3 is easily obtained by solving Eqs. Q 
and © to give /3e = which corresponds to a spherical shape. 
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The SO{6) limit 

For e = and X = Oj one recovers the the 50(6) limit 

H2 = -^Q(x)-Q(x) , X = . (11) 

Also in this case, the energy surface only depends on /3 

y2{P) - — (1+^2)2 • (12) 

The equilibrium value is given by /3e = 1 corresponding to a 7-unstable deformed shape. 
The SU{3) hmit 

For e = and x = -FV^ /2, the schematic Hamiltonian of Eq. reproduces the SU{3) limit 

H3 = -KQ{x)-Qix) , x = tIV7. (13) 

For this case, the potential-energy surface depends both on /? and 7 

kN{N - 1 



(l + /32)2 



Ap^±2^/2P^ COS37 + 



(14) 



The equilibrium values are given by Pe = and -fe = 0° for x = ~\"^ and by /3e = \/2 and 
7e = 60° for X = corresponding to prolate and oblate deformed shapes, respectively. 

2.4 Shape-phase transitions 

The transitional regions between any two of the dynamical symmetries of the IBM can be 
studied by combining the corresponding Hamiltonians, and varying the relative strength of the 
interaction terms. In 1980, Dieperink, Scholten and lachello showed that these transitional 
regions exhibit shape-phase transitions 

In 1996, the subject was revisited by Lopez-Moreno and Castahos j8| who investigated the 
shape-phase transitions for general (one- and two-body) IBM Hamiltonians using the full power 
of catastrophy theory TT. It was shown that the shape-phase diagram depends on two inde- 
pendent combinations of the parameters of the IBM Hamiltonian, called r\ and r2, which can 
be used to classify the equilibrium configurations. For the Hamiltonian of Eq. they are given 

by 



- + ^ - -^^V^ (15) 



In addition, closed analytic expressions were derived for the separatrix, which consists of the 
critical points 
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Table 1: Equilibrium values of the deformation parameters for the dynamical symmetry limits of the 
IBM. 



Limit e K X fi ^2 Pe 7e 



C/(5) 


e 







1 








SO{6) 





K 





-1 





±1 


SU{3) 





K 




4 

3 


|V2 


^2 0° 







K 




4 
3 


-|V2 


y/2 60° 



the spinodal and antispinodal points 

^ (9r| + 16)3/2 - 64 

ri. = -1 + ^ ' 2 , = . 17 

54r2 

In Tabled the values of ri and r2 are given for each of the dynamical symmetries discussed in 
the previous section. 

The methods used in ^ and jS| are closely related and give complementary results, which 
lead to a very simple and clear picture of the shape-phase diagram in nuclei as will be shown in 
the remainder of this section. 



The t/(5) — SO{Q) transitional region 

The transition between the spherical and 7-unstable deformed shapes can be studied by consid- 
ering the Hamiltonian 

H12 = H1 + H2 = ehd- KQ{<d)-Q{{)) . (18) 
The corresponding energy surface can be written as 

^^^^ Ne l + /?2 ^''(l+/?2)2' ^'^) 



where the energy functionals of Eqs. (|10j) and (|12j) were scaled by e and the number of bosons 
N. The transition is described by varying the control parameter rj = k{N — l)/e on the interval 
< 7/ < 00. The equilibrium value of /3 depends on -q 
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i.e. a spherical minimum for < 1/4 and a deformed one for rj > 1/4. The phase transition 
takes place at the critical point r/c = 1/4. Its nature can be determined by studying the ground 
state energy <?min(^) = ^iPe) and its derivatives at the critical point (Bj: if the ground state 
energy is discontinuous in the control parameter r] the transitition is zeroth order, and if the 
first (second) derivative is discontinuous it is first (second) order. The ground state energy of 
Eq. ((IH) is 

< ^ : '£'mm(??) = , 

1 (An - 1)^ 

V>-^ ■■ SmM = --^J^ ■ (21) 

The first derivative is continuous at r/c = 1/4, whereas the second derivative is discontinous. 
Therefore, the U{5)-SO{6) transitional region exhibits a second-order phase transition |^. 

An analysis in terms of catastrophy theory gives that the values of ri and r2 for the Hamil- 
tonian of Eq. ((THl) are given by 

1 - 4?7 , , 

= Fm? ■ = '''' 

In the r2-ri plane the transitional region for < ?7 < oo is represented by — 1 < ri < 1 and 
r2 = 0. The critical, antispinodal and spinodal points are obtained by taking the intersection 
with the curves corresponding to Eqs. (fT6|) and (|T7|) 

1 



r]c 



4 
1 



ric 


= 


= 


ria 


= r2a 


= 


ris 


= r2s 


= 



4 ' 

1 , , 

r]s = ris = r2s = Q . (23) 

This shows that the critical, antispinodal and spinodal points coincide which is a signature for a 
second-order phase transition, in agreement with the results obtained above from the analysis of 
the ground state energy 4 . In the left panel of Fig. ^ the potential-energy surface is given for 
three different values of the control parameter i]. For rj < rjc the equilibrium shape is spherical 
and for r] > i]^ it is deformed. There is no region of coexistence of spherical and deformed 
minima. 

The U{5) — SU{3) transitional region 

The transitional region between the spherical and axially deformed nuclei can be studied by 
considering the schematic Hamiltonian 

His = Hi+Hs = en, - kQ{x) ■ Q{x) , X = , (24) 

as a function of the control parameter r] = k{N — l)/e with < 77 < cxd. The corresponding 
energy surface is given by 



Vi{(5) + V3{Pn) r r, 



Ne I + /32 (l + /32)2 



4/32 ±2^2/3^ COS37 + 



(25) 
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Figure 1: Potential-energy surfaces, (left) The second-order C/(5) — 5*0(6) phase transition: rj < rjc 
(dashed), rj — r]c (solid) and rj > rjc (short-dashed), (right) The first-order U{5) — SU{3) phase transition: 
T] < rjs (dashed), i] ^ rjs (dotted), ij = rjc (solid), r] = i]a (dotted-dashed) and i] > rja (short-dashed). 



where the energy functionals of Eqs. (|lUjl and (|14j) are scaled by e and the number of bosons 
N. The analysis of the ground state energy as a function of r] is more complicated than in 
the previous case due to the presence of the (3^ cos 87 term. It was found that H13 exhibits a 
first-order phase transition at the critical point r/c = 2/9 (4j. 

An analysis with catastrophe shows that the values of the coefficients ri and r2 for H13 are 
given by 

from which it is easy to see that the transitional region between the U{5) and SU{3) limits is 
characterized by a straight line 

The critical, antispinodal and spinodal points are obtained by taking the intersection of the 
straight line of Eq. (|27|) with the curves corresponding to Eqs. and (|17j) . The critical point 
is characterized by 

2 1 8V2 , , 

rye = g , = > = • 

The equilibrium value of the deformation parameter at the critical point is given by |Hj 

i^c = = , (29) 



in agreement with For the antispinodal point one finds 



Va = ^ , ria = , = ±^y^ • (30) 
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The expressions for the spinodal point are rather comphcated, which is the reason why only 
numerical values are given 

7]s = 0.219 , ru = 0.075 , ra^ = ±0.748 . (31) 

In the right panel of Fig. ^ the potential-energy surface is given for five different values of the 
control parameter rj. For < r] < rjg the system only has a spherical minimum. At r] = r]s, it 
develops a second deformed minimum which becomes degerenate with the spherical minimum 
at T] = r]c- The spherical minimum disappears at r] = ija, and for r] > rja the system is deformed. 
For r]s < T] < rja, there is a coexistence of a spherical and a deformed minimum (see inset), 
the former is the lowest for r]s < r] < tJc, and the latter for r]c < V < Va- These features 
are characteristic of a first-order phase transitition. Recent experiments have provided strong 
evidence for ^^^Sm as a nucleus at the critical point of a phase transition between spherical and 
deformed shapes [Tl| . 



The 5*0(6) — SU{3) transitional region 

Finally, the transitional region between the SO{6) and SU{3) limits can be studied by 

H23 = H2 + H3 = -k'Q{0)-Q{0)-kQ{x)-Q{x) , x = ^\^ ■ (32) 

as a function of the control parameter ^ = k,/k' with < ^ < 00. The corresponding energy 
surface is given by 



N{N - 1)k' 
4/?2 



2^2 



(1 + /3 



2^2 



4/32 ± 2^/2/3^ cos 37 + ^/3^ 



(33) 



where the energy functionals of Eqs. H12|) and ()14|1 are scaled by A^(A^ — 1)k'. In this case, there 
is no phase transition |^. 

In terms of the coefficients ri and r2 one finds 



4+4e 



r2 



4eV2 



4 + 3C ' 4 + 3^ 

Also this transition corresponds to a straight line in the r2-ri plane 



-1± 



r2 



-- < n < -1 . 

3 - 



(34) 



(35) 



Since the line defined in Eq. ()35() does not intersect the curves for the critical, spinodal and 
antispinodal points, i?23 does not exhibit a phase transition. 
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Figure 2: Shape-phase diagram for the IBM as a function of ri and r2- The curves marked by c, s and 
a correspond to the critical, spinodal and antispinodal points, respectively. 

2.5 Shape-phase diagram 

In Fig. 121 the results of the previous section are summarized in terms of a shape-phase diagram 
of the IBM plotted as a function of the coefficients ri and r2 of Eq. ()15() . These coefficients 
are determined by the parameters of the Hamiltonian and completely classify the corresponding 
geometric shapes. A given Hamiltonian corresponds to a point in the r2-ri plane. However, 
note that the opposite is not true: a given point in the r2-ri plane corresponds not to a single 
Hamiltonian, but to a class of Hamiltonians which all share the same potential-energy surface, 
but differ in their kinetic energies. 

The transitional regions between the dynamical symmetries form a kite (or a double triangle) 
spanned by the straight lines derived for each one of them, see Eqs. ((221) i and (1^5]) . The 
dynamical symmetries occur at the end points of the kite, C/(5) (spherical), SO{Q) (7-unstable 
deformed) and SU{3) (prolate and oblate deformed). The critical points for the phase transition 
between spherical and deformded nuclei all lie on the curve defined by Eq. which is labeled 
by "c" in Fig. |21 The spinodal and antispinodal points are labeled by "s" and "a" , respectively. 

The region above the spinodal points (ri > ris) corresponds to spherical equilibrium shapes 
with only a spherical minimum. At the spinodal point (ri = rig), the system develops a second 
deformed minimum which becomes degerenate with the spherical minimum at the critical point 
(ri = Tie). At the antispinodal point (ri = ria), the spherical minimum disappears and the 
region below the antispinodal points (ri < ria) corresponds to deformed equilibrium shapes 
with only a deformed minimum. For rig < ri < ria, there is a coexistence of a spherical and a 
deformed minimum, the former is the lowest for rig < ri < ric and the latter for ric < ri < ria- 
These are the characteristic features of a first-order phase transitition. On the line r2 = 0, the 
spinodal, critical and antispinodal points coincide, corresponding to an isolated point with a 
second-order phase transition. The line r2 = and ri < has been identified with a first-order 
quantum phase transition for 7-soft nuclei [SI El • 

An often-used representation of the phase space of the IBM is the so-called symmetry or 
Casten triangle which was introduced in ^Hl as a schematic way to represent the properties of 
the Hamiltonian: the dynamical symmetries at the corners, the transitions between two of them 
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along the sides, and arbitrary combinations of the three of them in the inside of the triangle [2]. 
In Section [231 it was shown that the catastrophe theory treatment of [S] provides a derivation of 
this triangle, not as a schematic representation, but as a function of two control parameters ri 
and r2 which depend on the parameters of the Hamiltonian. This method is valid not only for the 
schematic Hamiltonian discussed in this contribution, but also for arbitrary IBM Hamiltonians. 



2.6 Shape-phase coexistence 

For first-order phase transitions between spherical and deformed shapes there is a very small 
region of coexistence. The right panel of Fig-^shows that the barrier that separates the spherical 
and deformed minima is very small, even in the inset it is hard to see. The height of the barrier is 
only a few percent of the excitation energy of the 2^ state. In a recent study, evidence has been 
presented for the existence of such coexisting phases at low energy in the transitional nucleus 

However, it is important to distinguish between this type of shape coexistence which occurs 
for spherical-deformed transitional nuclei and the far more common phenomenon of coexisting 
shapes due to different configurations I17j . The latter situations have been treated in an 
extension of the IBM in terms of a configuration mixing calculation, in which the normal and 
intruder configurations are described by IBM Hamiltonians and the configuration mixing by a 
coupling term ^H] 



^normal ^^ix 
^'mix -f^intruder 



(36) 



In this subsection, I briefly discuss some recent results obtained for the Ge isotopes [THl 1^. 
The even-even Ge isotopes show strong evidence for the coexistence of different geometric shapes. 
A sensitive probe for the configuration mixing is provided by two-neutron transfer reactions {p, t) 
and {t,p), in particular the ratio of cross sections for the excitation of the 0^ state and the ground 
state 

R = - ; . (37) 

In the absence of mixing, this ratio is zero {U{5) limit) or very small {SU{3) and S0{6) limits) 
[21 j . However, if the two configurations are strongly mixed, R can be different from zero. In 
Fig- El I show the results for {p, t) and {t,p) reactions between the Ge isotopes The enhance- 
ment in the cross sections around = 40 is an indictation of the importance of configuration 
mixing. 



2.7 Triaxial degrees of freedom 

The one- and two-body IBM-1 Hamiltonian can give rise to spherical, axially deformed and 
7-unstable deformed shapes. There is no stable triaxially deformed shape, unless one includes 
three-body interactions |22]. In the IBM-2 in which one distinguishes between protons and 
neutrons, there is in addition to the U{5), SU{3 and S0{6) limits, another dynamical symmetry 
called the SU*{3) limit, which has been associated with triaxial shapes {23 j . It is therefore 
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Figure 3: Comparison between experimental and theoretical cross sections for two-neutron transfer 
reactions {p, t) and {t,p) between Ge isotopes. N denotes the number of neutrons in the final Gc nucleus. 



of great interest to extend the discussion of phase transitions to IBM-2, since there are new 
transitional regions which are not present in IBM-1, namely U{5) — SU*{3), SU{3) — SU*{3) 
and 50(6) — SU*{3). Let's consider a schematic IBM-2 Hamiltonian 

H = e {Ud^ +nd,)- K [QjriXTr) + QuiXu)] ■ [QAXtt) + QuiXiy)] , (38) 

consisting of the d-boson energies and a quadrupole-quadrupole interaction. The SU*{3) limit 
arises for e = and Xvr = —Xu = /2. 

The analysis of the equilibrium shapes becomes a lot more complicated since in general one 
now has to study the properties of an energy surface that depends on 7 geometric variables: P-j^, 
l-K-, {f^u, lu) for the geometry of the proton (neutron) distribution and three Euler angles (G) = 
(Bi, 02, Qs) for the relative orientation of the proton and the neutron distributions. The general 
form of the potential-energy surface of IBM-2 can be found in Qi^. In the same reference it was 
shown that only IBM-2 Hamiltonians with a repulsive hexadecupole-hexadecupole interaction 
can give rise to so-called oblique shapes, i.e. an equilibrium shape with 7^ 0. Therefore, for 
the Hamiltonian of Eq. (|38p one can take = [251 126j . The energy surface is then given by 
the simplified form 



[2 

4cos(7^ - -fu) - 2y ^ XvPv cos(7,r + 27,.) 



"(l + /32)(l + /?2) 

[2 2 
-2y -j XttPtt cos{j^ + 27^) + - XtvXvPtvPu cos(27^ - 27j,) 



(39) 



where V{(3p,^p) is the energy surface for the proton {p = vr) and the neutron [p = v) part of 
the Hamiltonian, see Eq. @. 

The results that have been obtained so far for IBM-2 are mostly numerical |25|l26j: there is a 
second-order phase transition for U{h)-SU*{3) and SU{3)-SU*{3), whereas for SO{Q)-SU*{3) 
there is no phase transition. It would be very interesting to have analytic results for the phase 
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transitions between the four dynamical symmetries of the IBM-2. In this respect, the methods of 
catastrophe theory could be of help to obtain closed expressions for the critical and (anti)spinodal 
points, to provide insight into the order of the new phase transtions, and perhaps even to derive 
a shape-phase diagram for IBM-2, much in the same way as was done for the IBM-1 [S]. 

3 Random interactions 

A different method to study the geometric shapes of the IBM, and especially their robustness, 
is that of Hamiltonians with random interactions. 

Random matrix theory was developed to describe statistical properties of nuclear spectra, 
such as average distributions and fluctuations of peaks in neutron-capture experiments [271 128j . 
In this approach, the Hamiltonian matrix elements are chosen at random, while keeping some 
global symmetries, i.e. the matrix should be hermitean, and be invariant under time-reversal, 
rotations and reflections. Specific examples include the Gaussian Orthogonal Ensemble (GOE) 
of real-symmetric random Hamiltonian matrices in which the many-body interactions are un- 
correlated, and the Two-Body Random Ensemble (TBRE) in which the two-body interactions 
are taken from a distribution of random numbers j29j . For two-particle systems the two en- 
sembles are identical but for more than two particles, unlike in the case of GOE, in TBRE the 
many-body matrix elements are correlated. As a consequence, also the energy eigenvalues of 
states with different quantum numbers are strongly correlated, since they arise from the same 
Hamiltonian. 

The latter aspect was investigated in shell model calculations for even-even nuclei in the sd 
shell and the pf shell [H^. An analysis of the energy spectra of an ensemble of random two- 
body Hamiltonians showed a remarkable statistical preference for ground states with angular 
momentum and parity = 0'*', despite the random nature of the two-body matrix elements 
both in sign and relative magnitude. A similar preponderance of 0"*" ground states was found 
in an analysis of the IBM with random interactions [HJ. In addition, in the IBM evidence 
was found for both vibrational and rotational band structures. According to the conventional 
ideas in the field, the occurrence of L = ground states and the existence of vibrational and 
rotational bands are due to very specific forms of the interactions. Therefore, these unexpected 
and surprising results with random interactions have sparked a large number of investigations 
to explain and further explore the properties of random nuclei (see e.g. the reviews by Bijker 
and Frank 32 , Zelevinsky and Volya j33j and Zhao, Arima and Yoshinaga |34j). 

In this section, I discuss the phenomenon of emerging regular spectral features from the 
IBM with random interactions and its relation with the underlying geometric shapes and critical 
points. 

3.1 Regular features and robust results 

In order to study the geometric shapes associated with the IBM, let's consider the most general 
one- and two-body IBM Hamiltonian 

(40) 



H 



1 



1 



N - 1 



-Ho 
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where Hi^b contains the boson energies 



2 



2^/2 
1 



Hi^b = ess^ -s + ead^ -d , (41) 
and H2^b the two-body interactions 

+ E CA^((i^xdt)(A).(Jxd)(^) 
A=0,2,4 

1= [(dt X dt)(o) . {s X + (st X st)W . (J X dfy 
'{d^ X . (d X + (st X dt){2) . (J X . (42) 

The nine parameters of the IBM Hamiltonian of Eqs. (|4U|) - (|42|) . altogether denoted by 

(x) = {es,ed,uo,U2,co,C2,C4„vo,V2) , (43) 
are taken as independent random numbers Xi {i = 1, . . . ,9) on a Gaussian distribution 

P{x,) = e-^'?/2-'/a^/2^ , (44) 

with zero mean and width a. In this way, the interaction terms are arbitrary and equahy hkely 
to be attractive or repulsive. The spectral properties of each Hamiltonian are then analyzed by 
exact numerical diagonalization |31j . 

The properties of the energy spectra can be studied by the energy ratio 

In Table m the values of R4/2 are presented for each of the dynamical symmetries: R4/2 = 2.0, 
2.5 and 3.3 for the U{5), SO{6) and SU{3) limits, respectively. Fig. |3] shows a remarkable 
result: the probability distribution P(i?4/2) has two very pronounced peaks, one at ~ 1.95 and 
a narrower one at ~ 3.35. These values correspond almost exactly to the harmonic vibrator and 
rotor values (see the results for the U{5) and SU{3) limits in Table No such peak is observed 
for the 7-unstable or deformed oscillator case (50(6) limit). This can be understood from the 
fact that 7 dependence of the potential-energy surface arises from a single parameter, V2- For 
a 7-unstable solution to occur, V2 has to be zero. For V2 ^ the energy surface has a prolate 
or an oblate minimum (whether the global minimum is spherical or deformed depends on the 
other terms in the Hamiltonian as well). 

Energies by themselves are not sufficient to decide whether or not there exists band structure. 
Levels belonging to a collective band are connected by strong electromagnetic transitions. In 
Fig. a correlation plot is shown between the ratio of B{E2) values for the 4^ 2f and 
2f of transitions and the energy ratio R4/2- Table I^J shows that in the large limit this 
ratio of B{E2) values is 2 for the harmonic oscillator (C/(5) limit) and 10/7 for the deformed 
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Table 2: Energies of B{E2) values in the dynamical symmetry limits of the IBM for the Hamiltonians 
of Eqs. ®, CU and (O. 





E{4+)-EiO+) 
E{2+)-E{0+) 


EiO+)-E(0+) 
E{2+)-E{0+) 


B{E2;A+~*2+) 
B{E2;2+'^0+) 






2 


2 


2{N-1) 

N ~^ 


2 


50(6) 


5 
2 


+ 1 


10{Af-l)(iV+5) 
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oscillator and the rotor {S0{6) and SU{3) limit, respectively). The right panel of Fig. |2 shows 
a strong correlation between the first peak in the energy ratio and the vibrator value for the 
ratio of B{E2) values (the concentration of points in this region corresponds to about 50 % of 
all cases), as well as for the second peak and the rotor value (about 25 % of all cases). 

Despite the random nature of the interaction strengths both in relative size and sign, the 
ground state still has L = in the majority of the cases. Fig. Elshows the percentages of ground 
states with L = and L = 2 as a function of the boson number N (solid line). One sees a clear 
dominance of ground states with L = with ~ 60-75%. For N = 3k (a. multiple of 3) one sees 
an enhancement for L = and a decrease for L = 2. The sum of the two hardly depends on the 
number of bosons. 




Figure 4: Probability distribution P(i?4/2) of the energy ratio -R4/2 (left) and correlation between ratios 
of B{E2) values and energies (right) for the IBM with random one- and two-body interactions. The 
number of bosons is iV = 16. 



14 



100 



L=0 



80 



M 60 




^ 40 - 



L=2 



20 








6 



8 



10 



12 



14 



16 



N 



Figure 5: Percentages of ground states with L — (top) and L — 2 (bottom) in the IBM with random 
one- and two-body interactions as a function of the number of bosons N: calculated exactly for 10,000 
runs (solid line) and in mean- field approximation (dashed line). 

3.2 Mean-field analysis 

The results obtained in studies of the nuclear shell model and the IBM with random interactions 
are very surprising and unexpected in the sense that, according to the conventional ideas in the 
field, the occurrence of L = ground states and the existence of vibrational and rotational 
bands are due to very specific forms of the interactions. 

The basic ingredients of the numerical simulations are the structure of the model space, the 
ensemble of random Hamiltonians, the order of the interactions (one- and two-body), and the 
global symmetries, i.e. time-reversal, hermiticity and rotation and reflection symmetry. The 
latter three symmetries cannot be modified, since nuclear levels have real energies and good 
angular momentum and parity. It has been shown that the observed spectral order is a robust 
property that does not depend on the specific choice of the ensemble of random interactions 
pt). 35, 36) , the time-reversal symmetry |35j , or the restriction of the Hamiltonian to one- and 
two-body interactions ; 37^ . These results suggest that that an explanation of the origin of the 
observed regular features has to be sought in the many-body dynamics of the model space and/or 
the general statistical properties of random interactions. 

In this section, the origin of the regular features that emerge from random interactions in 
the IBM |38j are studied using standard Hartree-Bose mean-field methods in which the trial 
wave function is the boson condensate of Eqs. and © (see Section [221) • 

Vibrations 

The potential-energy surface associated with the IBM Hamiltonian of Eqs. H40() - H42() is given by 
its expectation value in the coherent state 




(3' 
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a3/3^ cos 37 -|- 04/3'^ 



(46) 
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where the coefficients Oj are hnear combinations of the parameters of the Hamiltonian 
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(47) 



The energy surface of Eq. ()46|) provides information on the distribution of shapes that the 
model can acquire. The values of /5e and 7e that characterize the equilibrium configuration of 
the potential-energy surface depend on the coefficients 04, and 02- The equilibrium shapes 
can be divided into three different classes: 

• f3e = 0: an s-boson or spherical condensate 

• < /5e < 00 with 7 = 0° or 7 = 60°: a deformed condensate with prolate or oblate 
symmetry, respectively, and 

• Pe = 00: a d-boson condensate. 



Each equilibrium configuration has its own characteristic angular momentum content. Even 
though the angular momentum states are not projected from the coherent state, the ground 
state angular momentum can be obtained from the rotational structure of the condensate in 
combination with the Thouless-Valatin formula for the corresponding moments of inertia. This 
procedure is described in detail in j^H]. The application to the IBM with random interactions 
leads to the following results. 

• The s-boson condensate corresponds to a spherical shape. Whenever such a condensate 
occurs (in 39.4 % of the cases), the ground state has L = 0. 

• The deformed condensate corresponds to an axially symmetric deformed rotor. The or- 
dering of the rotational energy levels L = 0, 2, ... , 2N is determined by the sign of the 
moment of inertia 



The moment of inertia T3 depends in a complicated way on the parameters in the Hamil- 
tonian. The deformed condensate occurs in 36.8 % of the cases. For X3 > the ground 
state has L = (23.7 %), while for X3 < the ground state has the maximum value of the 
angular momentum L = 2N (13.1 %). 



Rotations 



'rot — 



(48) 
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• The d-boson condensate corresponds to a quadrupole oscillator with quanta. Its rota- 
tional structure has a more complicated structure, since the rotational excitation energies 
depend on two moments of inertia 



which are associated with the spontaneously broken three- and five-dimensional rotational 
symmetries of the d-boson condensate. The d-boson condensate occurs in the remaining 
23.8 % of the cases. The results in Table 01 can be understood qualitatively as follows. For 
T5 > the ground state has r = for N even or r = 1 for N odd (~ 4 %), while for 
X5 < the ground state has the maximum value of the boson seniority t = N (~ 19 %). 
For r = and r = 1 there is a single angular momentum state with L = and L = 2, 
respectively. For the t = N multiplet, the angular momentum of the ground state depends 
on the sign of the moment of inertia X3 . For X3 > the ground state has L = for N = 3k 
or L = 2 for N ^ 3k (9 %), while for T3 < the ground state has the maximum value of 
the angular momentum L = 2N (10 %). 

Table 01 shows that the spherical and deformed condensates contribute constant amounts of 
39.4 % and 23.7 %, respectively, to the L = ground state percentage, whereas the contribution 
from the d-boson condensate depends on the number of bosons N. The L = 2 ground states 
arise completely from the d-boson condensate solution. Fig. shows the percentages of ground 
states with L = and L = 2 as a function of the total number of bosons A^. A comparison of 
the results of the mean-field analysis (dashed lines) and the exact ones (solid lines) shows an 
excellent agreement. There is a dominance of ground states with L = for ~ 63-77 % of the 
cases. Both for L = and L = 2 there are large oscillations with which are entirely due to 
the contribution of the d-boson condensate. For N = 3k one sees an enhancement for L = 

Table 3: Percentages of ground states with L = 0, 2 and 2A^, obtained in a mean-field analysis of the 




(49) 



IBM. 



Shape 
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and a corresponding decrease for L = 2. In the mean- field analysis, the sum of the two which 
accounts for ~ 77 % of the cases, hardly depends on the number of bosons, in agreement with 
the exact results. For the remaining ~ 23 % of the cases, the ground state has the maximum 
value of the angular momentum L = 2N. 

In conclusion, the emergence of regular features from the IBM with random interactions 
can be explained in a Hartree-Bose mean-field analysis of the random ensemble of Hamiltoni- 
ans, in which different regions of the parameter space are associated with particular intrinsic 
states, which in turn correspond to definite geometric shapes |38| . There are three solutions: a 
spherical shape carried by a single state with L = 0, a deformed shape which corresponds to 
a rotational band with L = 0, 2, ... , 2N, and a condensate of quadrupole bosons which has a 
more complicated angular momentum content. The ordering of rotational energy levels depends 
on the sign of the corresponding moments of inertia. The mean-field analysis explains both the 
distribution of ground state angular momenta and the occurrence of vibrational and rotational 
bands. The same conclusions hold for the vibron model for which a large part of the results has 
been obtained analytically [10] . 

3.3 Evolution of quasi-beta energies 

Recently, an interesting empirical correlation was pointed out [IJ between of Eq. (|l5|l and 
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The left panel of Fig. |H1 shows the experimental values for the nuclei with neutron number 
84 < < 96 and the right panel the results of a IBM calculation in which the parameters of 
the CQF Hamiltonian of Eq. ^ are taken randomly, but are restricted to the physically allowed 
region, i.e. e > 0, k > and — \/7/2 < x ^ V7/2. For vibrational nuclei both energy ratios 
are equal to 2, since the O2 ^1 states belong to the two-phonon multiplet, whereas for 
rotational nuclei the 4^ state is a rotational excitation and the state a vibrational excitation 
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Figure 6: Correlation plot between the energy ratios R0/2 and -R4/2: (left) experimental values HT" and 
(right) theoretical values calculated in the consistent Q-formulation of the IBM. 
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Table 4: Comparison of energy ratios for the critical points of the IBM Hamiltonians H12 and Hi^ with 
= 11 and the critical point symmetries 



C/(5) - SO{<o) 
Vc = 1/4 



U{5) - SU{3) 
ric = 2/9 



X(5) 




2.22 
2.61 



2.20 
3.03 



2.33 
2.41 



2.91 
5.67 



(see Table 12) • Fig. |H1 shows that this trend is satisfied both by the data and by schematic IBM 
calculations in which only some minimal restrictions are imposed upon the Hamiltonian. 

3.4 Signatures of critical point nuclei 

Finally, I briefly address the question of experimental signatures of critical point nuclei |10j . 
Phase transitions in nuclei can be tested experimentally by measuring observables that are 
particularly sensitive to them, such as separation energies 8211-, isomer shifts (5(r^) and the 
transition rates B{E2; 2^^ —>■ O2) Other observables that have been suggested as measures 
of critical point nuclei are the energy ratios R4/2 and i?o/2- Table |1] shows the values of these 
ratios for the IBM Hamiltonians H12 and H13 at the critical points r/c = 1/4 and r/c = 2/9, 
respectively. A comparison with the values for the critical point symmetries, E{5) and X(5) 0, 
shows that the energy ratios for the critical points in the IBM are rather different from the those 
for the critical point symmetries. This is not so surprising, since the potential-energy surfaces 
for IBM Hamiltonians at the critical points are not square-well potentials (see Fig. ^ as in the 
critical point symmetries [2]. Moreover, one has to take into account that there is an entire class 
of IBM Hamiltonians which all share the same potential-energy surface and hence have the same 
critical points, but differ in their kinetic energies. Therefore, there is not a single value of the 
energy ratios R4/2 and -R0/2 that can be used to identify a nucleus as a "critical point" nucleus. 
Instead one should rather study the change in energy ratios, B{E2) values and/or quadrupole 
moments in a chain of nuclei to determine whether these nuclei exhibit a shape-phase transition. 

4 Summary, conclusions and outlook 

In these lecture notes I have presented a review of shape-phase transitions in nuclei and random 
interactions. Maybe at first sight, these two topics have little to do with one another, but both 
were shown to provide evidence for the existence of robust shapes in collective nuclei. The 
transitions between different shapes occur very rapidly, typically with the addition or removal 
of one or two pairs of nucleons. 

An analysis in IBM-1 and IBM-2 showed that the phase transtions between spherical and 
axially deformed nuclei are of first order, whereas those between spherical and triaxially deformed 
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and between axially deformed and triaxially deformed are of second order. 

Two-neutron transfer reactions {p,t) and {t,p) for the Ge isotopes were shown to provide 
a sensitive probe of shape-phase coexistence of a normal and an intruder configuration. In 
particular, the ratio of cross sections for the excitation of the first excited 0+ state and the 
ground state is either very small or zero in the absence of configuration mixing, but may be 
enhanced for strongly mixed configurations. 

The importance of the concept of geometric shapes in nuclei was confirmed by a study of the 
IBM with random interactions. Despite the random nature of the interactions a highly surprising 
degree of regularity was observed, in particular the dominance of ground states with L = and 
the occurrence of vibrational and rotational band structures. It was shown that the geometric 
shapes associated with IBM Hamiltonians play a crucial role in understanding the origin of 
these regular properties. These shapes are a reflection of an intrinsic geometry associated to 
the many-body dynamics of the model space which is sampled by the statistical nature of the 
random interactions, but which is quite independent of them. The observed spectral order is a 
robust property that arises from a much larger class of Hamiltonians than is usually thought. 
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